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Abstract 

We study pairing correlations in ultrasmall superconductor in the nanoscopic 
limit by means of a toy model where electrons are confined in a single, multiply 
degenerate energy level. We solve the model exactly to investigate the tem- 
perature and magnetic field dependence of number parity effect (dependence 
of ground state energy on evenness or oddness of the number of electrons). 
We find a different parity effect parameter to critical temperature ratio (~4 
rather than 3.5) which turns out to be consistent with exact solution of the 
BCS gap equation for our model. This suggest the equivalence between the 
parity effect parameter and the superconducting gap. We also find that mag- 
netic field is suppressed as temperature increases. 
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Single electron tunnelling experiments of Ralph, Black and Tinkham§ (RBT) initiated 
interest on the possibility of superconductivity in a nanometer size Al grain. Although ex- 
istence of superconductivity in a very small metallic particle has been discussed long time 
ago! it has been possible to investigate such small systems only recently with experiments of 
RBT. Anderson's criterion, that superconductivity ceases when discrete energy level spac- 
ing becomes comparable with superconducting gap is correct in general, should be treated 
carefully at such small scales. 

Instead of macroscopic properties like zero resistance and Meissner effect, superconduc- 
tivity in nanometer size metallic grain manifests itself as an odd-even or number parity effect, 
i.e. dependence of physical properties on whether grain has odd or even number of electrons 
in it. Evaluation of discrete energy spectrum using BCS model or exact diagonalization of 
finite systems has been the main subject of theoretical studiesihi. (For a review see Ref.i). 

For parabolic energy dispersion, which is the case for a grain with perfect symmetrv. 
discrete energy levels are multiply degeneratei^. Figure 25.1 of Ref.i and Fig. 1 of Ref.0 
show explicitly that for parabolic dispersion degeneracy of energy levels can be very large. 
For example, the number of solutions satisfying = nf + + n"^, all being integer, can 
approach to 40 for n < 100. Although disorder and irregular shape split the degeneracy, 
as long as they are weak, we can still talk about near degeneracy. The main point is that 
typical scale of dimensional energy quantization 61 = hvp/L is, L being the size of the 
grain, much larger than level spacing 62 — Ep/N, where Ep is fermi energy and is 
the total number of electrons, even for disordered samples. If attractive electron-electron 
interaction is much less than di, but yet larger than 62, we can approximate the system 
as a single, multiply degenerate energy level provided that disorder is not too strong. A 
single energy level is certainly drastic approximation. Our aim is to find out possibility 
of superconductivity, or pairing correlations, in this extreme limit. Physical systems are 
probably between bulk BCS model and the nanoscopic (or quantum) limit. By examining 
the behavior of a superconductor at both ends we expect to figure out properties of real 
systems. Finally, our earlier worksi'0 show that exact solution of multi-level model agrees 
very well with single energy level model in the large level spacing (or nanoscopic) limit. 
The formula, obtained in the single level approximation, for the ratio of the number parity 
effect parameters for odd and even cases gives results within one percent neighborhood of 
the exact solution when level spacing is about ten times larger than the superconducting 
energy gap. 

Our starting point is a single, multiply degenerate energy level. We consider d states 
and assume that electrons in the same state are paired. Therefore the model Hamiltonian 
becomes 

d 

H = -g 4'T4'iO«i««T7 (1) 

n,n'=l 

where aj^^ (ono-) creates (annihilates) an electron in state n with spin cr. Such a model has 
already been used long time ago by Mottelson to study pairing of nucleons^. In our case 
n describes a state where time-reversed partner is itself. Hence, we assume that electrons 
occupying the same state are paired and they can be scattered into anyone of the unoccupied 
states together with the same amplitude. It is an immediate consequence of the form of the 
interaction that singly occupied states give vanishing matrix elements. Single electrons 
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simply block a state from any scattering event. Therefore, their mere effect is to eliminate 
certain states in the Hamiltonian. On the other hand, doubly occupied and empty states can 
be denoted in pseudo-spin representational which allows us to diagonalize the Hamiltonian 
easily. 

Pseudo-spins s are introduced by the definition s„+ = a^-^al^ and Snz = ^{O'n^O'ni + 
al^ttni — 1) which allows us to write the Hamiltonian as if = —gS-^-S- where S = J2n is 
the total spin and 5*^ are corresponding raising and lowering operators. Finally, when it is 
written as if = —(7(8^ — 5*^ + Sz), we obtain energy eigenvalues as 

E{s) = -j{N - s){2d- s - N + 2) (2) 



where 



0,2,4, ■■■,2rf-A^,---,A^ for N even , . 

1, 3,5, 2d- N,---,N for N odd. 



Here, seniority number s = d — 2S has been introduced to simplify the expressions. We 
note that for a given number of electrons N, z-component Sz of the pseudo-spin is fixed and 
hence eigenvalues are labelled by only one quantum number S. Degeneracy of each S, i.e. 
several different ways of obtaining the same total spin value, makes the problem nontrivial. 
Therefore, we have have to find the total number of states Q{s) having the same seniority 
number s to evaluate thermodynamical quantities. Let us start with the ground state which 
corresponds to the largest possible total spin 5* or equivalently the smallest possible seniority 
s. For even A^, s = state is unique whereas for odd N, due to un-paired electron there is 
2d-io\d degeneracy. Here, we assume that the grain interacts with a heat bath where spin 
flip is possible which gives the factor 2. 

Next, we consider the excited states which can be created not only by pair breaking 
(single particle excitations) but also by changing total spin value with fixed number of 
pairs (collective excitations). Single particle and collective excitations turn out to have the 
same energy spectrum. Hence, the total degeneracy Q{s) can be obtained by summation 
of degeneracies of each configuration. Contribution of un-paired electrons to degeneracy is 
nothing but the number of combinations of d states taken s at a time. Including spin degree 
of freedom results in a factor 2*. Different ways of obtaining the same total spin S with 
spin-1/2 particles, which is given by 

(4) 

brings an additional contribution. Total degeneracy is found by proper use of the two 
formulae. Summation over all possible configurations having the same seniority gives us the 
total degeneracy 

(5) 








where I is equal to |, -^^-^ for even and odd A^, respectively0. 
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Before going to evaluate thermodynamical quantities let us calculate the ground state 
energy as a function of number of electrons A^, for a given degeneracy d. Figure. 1 shows 
d =20 case. For N =20, the level is half-filled and we reach the minimum. This is an 
expected result since both the number of pairs and the number of the empty states into 
which they scatter are maximized to have the optimum energy. Zig-zag structure of the 
curve is due to the parity effect. When is odd there is a unpaired electron which is 
neither scattered out and nor allows any pair to be scattered into the state it occupies. This 
blocking effect effectively decreases the number of states by one. 

Noting the fact that Zeeman energy shifts of discrete levels are much more important than 
orbital diamagnetismlll we can immediately write down energy eigenvalues in the presence 
of a uniform magnetic field as 

Eh = Eh=o - fiBH{n^ - rii), (6) 

where and denote the number of spin-up and spin-down electrons, respectively. Al- 
though, we find Eh easily, degeneracy expression is no more simple but instead changes with 
magnetic field. However, for small number of states it is not difficult to count the energy 
levels to find degeneracy of each energy eigenvalue. 

Having degeneracy of each state we can evaluate specific heat Cy = as a function of 
temperature where 

Observation of specific heat variation is one of the most direct ways to find critical temper- 
ature Tc and magnetic field He- 

Figure. 2 shows variation of Cy{T) with temperature for different number of electrons. 
Discontinuity in the specific heat, which is a indication of a second order phase transition 
in macroscopic systems, now manifests itself as a relatively fast variation in Cy{T) curve. 
Absence of a sharp change results in ambiguity in determination of T^. We are going to 
identify the temperature at which the maximum of Cy{T) curve is reached as T^. We note 
that the linear background, i.e. contribution to Cy{T) which is of the form const. T, is 
absent. This is because linear term is a result of Fermi-Dirac distribution and continuous 
spectrum near Fermi energy. Since in our case, we have a single level, Cy{T) curve has a 
flat background. 

In Fig. 3 we plot critical temperature for different magnetic fields as a function of number 
of electrons. Higher value of Tc for even number of electrons is an indication of number parity 
effect which can be defined as 



(8) 

where E"^^ is the ground state energy for 2A^ electrons. Ap has no label because for our 
model it can be shown that Ap = gd/2, independent of A^. However, at finite temperature 
Ap should be redefined. One possibility is to replace ground state energies by thermal 
averages where we obtain Ap(T) as a decreasing function of temperatures. We note that 
for H = 0, kTc/g ^ 2 while Ap/g = 5 which gives 2Ap/kTc = 5. With increasing d we 
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observe that ratio decreases. For d = 100, which is not plotted due to invisibihty small 
odd-even dependence, 2Ap/kTc — 4.6. Here the interesting point is that BCS gap equation 



can be solved analytically, with V]^^ = —g, E-^^ = E and Aj^. = A, to give kTc/g = d/4 ~ 
Ao/2g(, and we find 2Ao/kTc = 4 (rather than 3.5). 

Finally, in Fig. 4, we plot critical magnetic as a function of temperature. We attribute 
the feature at low temperatures to slightly ambiguity in the definition of critical temperature. 
In interpreting Hc{T) curve we should keep in the mind that the critical magnetic field is 
determined by Zeeman splitting. For a different geometry, for example for a ring shaped 
conductor, orbital coupling can also play an important role. 

In conclusion, we can understand basic physics of thermal and magnetic properties of a 
nanoscopic superconductor using our very simple model which can be solved exactly. In the 
extreme nanoscopic limit where a single degenerate or nearly degenerate energy level is left 
we predict a different gap to critical temperature ratio. Our results, with increasing number 
of levels, approach to exact solution of the BCS gap equation which gives 2Ao/A;Tc = 4 
(rather than 3.5). This suggest the equivalence between the parity effect parameter Ap and 
the superconducting gap A. Finally, we observe that the critical magnetic field is suppressed 
with increasing temperature. 
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FIGURES 



FIG. 1. Ground state energy with respect to the number of particles in a grain with d = 20. 

FIG. 2. Specific heat of a grain with d = 100 and N = 100, 50, 2. Smooth transition occurs 
from pair-correlated region to normal region around the maximum of Cy. 

FIG. 3. Critical temperature with respect to number of particles for a grain with d = 10 for 
different magnetic fields ixH/g = 0,2,3. Number parity dependent Tc decreases with increasing 
magnetic field. 

FIG. 4. Critical magnetic field with respect to temperature where d = 10 and = 10. He is 
suppressed with temperature. 



6 



300 



200 - 



100 - 












10 



20 



30 



kT/g 



